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Abstract 

The description of a three-dimensional Ising-hke magnet in the presence 
of an external field in the vicinity of the critical point by the collective vari- 
ables method is proposed. Using the renormalization group transformations, 
the scaling region size is defined as a function of temperature and field. The 
obtained expressions for the free energy, equation of state and susceptibil- 
ity allow one to analyse their dependence on microscopic parameters of the 
system. The critical exponents of the correlation length and order param- 
eter are calculated as well. The results agree qualitatively with ones ob- 
tained within the framework of the parametric representation of the equation 
of state and Monte-Carlo simulations. The calculations do not involve any 
parametrization, phenomenological assumptions and adjustable parameters. 
The approach can be extended to models with a multicomponent order pa- 
rameter. 



1 Introduction 

The description of the phase transitions in the three-dimensional {3D) magnets, 
usually, is associated with the absence of exact solutions and with many approximate 
approaches for obtaining different system characteristics Ij. Even the simplest Ising 
model, despite the variety of works devoted to its study, is still of interest [2]. 
Inclusion of the external factors affecting on the system behaviour complicates the 
problem. 

Some of the more popular approaches, used for the study of critical phenomena, 
are based on the series expansions of different characteristics in small parameters, 
scaling theory, renormalization group (RG) transformations, etc. 

Let us recall two kinds of the RG theory. One of them is based on the pertur- 
bation theory (see, for example, 0). The non-Gaussian distributions are used in 
another one Q . An advantage of the non-perturbative methods is the possibility to 
avoid so-called infrared (IR) problems, which exist in the pertubative approaches. 



One of the most popular non-perturbative approaches employs the exact RG equa- 
tion (see jUiniEl)- The application of the non-perturbative methods is effective for 
investigating the Ising universal behaviour in high energy physics. In |7] , the critical 
exponents for systems in the universality class of the 3D Ising model are computed 
in the local potential approximation (LPA), that is, within the framework of the 
Wegner-Houghton equation. The similar calculations are performed for different IR 
cutoffs, which can be chosen freely jH]. In both papers [7] and [Hj, the effective 
potential contains the odd part, which allows one to calculate the critical exponents 
related to the field dependence. 

We propose the description of the critical behaviour of a 3-D Ising-like magnet 
in the presence of an external field by some alternative non-perturbative approach 
called collective variables (CV) method OITO]. It exhibits some similar features 
of the Wilson exact RG approach. In particular, the Kadanoff hypothesis for the 
construction of the effective lattice, which is commensurate with the size of long- 
range interactions, is used in both methods. In contrast to other methods, the 
calculations carried out within the CV approach are performed on the microscopic 
level without any phenomenological assumptions and adjustable parameters. Star- 
ting from the Hamiltonian, the CV method makes it possible to obtain the critical 
exponents as well as the thermodynamic and structural system characteristics. In 
particular, we get the free energy, equation of state and susceptibility as functions 
of temperature and field. 

Using this method, we have already obtained thermodynamic characteristics in 
the form of a power series in the scaling variable in the regions of the weak and strong 
fields for temperatures above and below Tc (Tc is the phase transition temperature 
in the absence of an external field) [TTl I12j . In the mentioned papers as well as in 
jTB] , the size of the scaling region has been defined only by the temperature or by the 
field. In additional to even powers of CV, the effective Hamiltonian has contained 
the odd part represented by the linear and cubic terms, which are related to the 
presence of the field. In contrast to the works |llnl2j, the present calculations are 
carried out without using power series in scaling variable. The size of the scaling 
region is a function of both the temperature and field variables. Due to removing 
the cubic term (by the corresponding substitution of variables), the odd part of the 
effective Hamiltonian is represented only by the linear term. The obtained results 
are valid in the whole field-temperature plane. 



2 RG transformations 

The system Hamiltonian 

^ = -^E*(nj)^i^j-/iE^i (1) 

includes the interaction potential $(rij) = Aexp (— ry/fc), which is the function of 
the distance between spins (Ji located at sites 1 and j of a simple cubic lattice. The 
microscopic parameter b is the radius of effective interaction, A is a constant. The 



second term on the right-hand side of formula (0) is related to the presence of an 
external magnetic field h. 

The CV pk are introduced by means of the functional representation for operators 

Pk = (v^)-^E^iexp(-^kl). (2) 
1 

Using the Jacobian of transition from the set of N spin variables a\ to the set of CV 
Pk and the approximation of the Fourier transform of the interaction potential 

^^''^ \ $0 = $(0)$, Bo<k<B, ^''^ 
we start from the partition function in the following form: 

(dp)^" exp[-aiVNoPo - - ^ d{k)p]^p_k 

^ k<Bo 

Pki--Pk4^ki+...+k4]- (4) 



Here $ is the small constant, (5ki+...+k4 is the Kronecker symbol, Nq = Nsq , Bq = 
B/sq. The quantity B = ir/c determines the first Brillouin half-zone, c is the 
lattice constant, and the initial division parameter sq defines the region, where the 
parabolic approximation for $(fc) is valid [see Eq. The CV pk are associated 

with modes of spin moment density oscillations. The variable po corresponds to 
the system order parameter. In the general case, the expression in the exponent of 
relation (jU is the n-power polynomial, where n takes on the integer values. The 
polynomial with coefficients a„ determines Jacobian of transition to the set of CV 
iniCIIj. Performing our calculations, we are restricted by the simplest non-Gaussian 
distribution of the order parameter fluctuations with n = 4 (p'^ model). In order 
to cancel the cubic term, the corresponding substitution of variables is carried out. 
The CV method allows us to perform a well controlled calculation at each order of 
the truncation. In this respect, the term "model" may be associated with the order 
(models p^, p^, etc.). 

The initial coefficients in (@]) are expressed by relations 

a^ = -s^J^h', d(fc) =a2 + /3$o-/3$(A;), 

a2 = 1 - 3c;j, 04 = 6Cf^. (5) 

Here h' = (3h is the dimensionless field, (3 = 1/kT is the inverse temperature {k is 
the Boltzmann constant), and c)^ = 50^^(1 — 4/i'^)/3. The first term in the exponent 
of the partition function (0)) is related to the presence of an external field. The 
coefficients 02 and 04 depend on the field weakly. 

The calculation of the partition function (j3]) is performed using the step-by-step 
integration with respect to the CV. The procedure of such an integration is described 
in detail in HDl El- For this purpose, we divide the phase space of the CV pk 
into layers with the division parameter s. The integration begins from the variables 



with large values of the wave vector and terminates at pk with k 0. In each 
nth layer, which corresponds to the region of wave vectors -Bn+i < k < B^, where 
Bn+i = Bn/s, the Fourier transform of the potential is replaced by its average 
value. For simplicity of estimations, we do not take into account the correction 
for the potential averaging. Such an approximation corresponds to the LPA in the 
Wilson exact RG approach. Inclusion of this correction leads to an appearance of 
the critical exponent t], which characterizes the behaviour of the pair correlation 
function for T = Tc. Such calculations have been performed within CV method 
in the absence of an external field. The formal part of the procedure has already 
been presented in detail for the model in j^l . The quantity rj takes on the 
value 7] 0.024 For comparison, the exponent rj = 0.044 was obtained using 
the quartic effective potential and method based on an exact flow equation for a 
coarse-grained free energy Including the above-mentioned correction and the 
related shift of the fixed point does not qualitatively change the main thermodynamic 
characteristics of the system. 

Thus, as a result of performing n iterations, the partition function is presented 
in the form 

Z OC QoQl---QnIn+l, (6) 

where quantities 

Qn = [QiP^''-'^)Qidn)f- (7) 

are the partial partition functions of the nth block structure. Here 

Qidn)=(^]\l{l-W2^) (8) 



as well as 



pr"-(^)%(i+«4'-') 

(n) _ /^drJyBn+l,Bn) 



hr = V6 "t " . (9) 



The coefficients 7, 71, Gq, ^2 and G2 are independent of the field and are given in 
The quantity dn{Bn+i, Bn) = dn{0) + qs~'^'^ is determined by the averaged Fourier 
transform of the interaction potential Bn) with the help of the relation 

dn{k) = + mSn+l, Bn) - mk). (10) 

Here q = qP^{0) and q = (bn / c)"^ s^'^ {1 + The quantity In+i in Eq. ® has the 
following form: 



|(rfp)^"+^exp|-aS'^+^W.+iPo 



^ C?n+l(A;)PkP-k 

k<B„+i 
Ari+1) 



X! Pki---Pk4'^ki+...+k4 ^ 



(11) 



The coefficients af^^^^ are defined by recurrence relations (RR) 
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4"^'^ = /oo(4"V/'(i + «2/^^"^ 

at^'^ = s-''U,at\l + a,h^'^), (12) 

where quantities /oo, 0-2^ foi and 04 are constants JA^. Performing the scaling 
transformations pk = -sp^ with introduction of the notations Un = s^ai^\ r„ = 
s'^"'dn{0), Un = s'^^'oi'' and taking into account the expression for /i2"'* from Eqs. 0, 
the RR ()12|) can be presented in the matrix form 

Tn-r* I . (13) 

Un - U* 

Here Rn = s'^+^Z^ R22 = sVooa2A/6, i?24 = sVoo^^*"^''^ ^42 = 5/010476^*^/^ and 
-R44 = s/oi- The coordinates of the fixed point w* = 0, r* = —q and u* = q{l — 
■s~^)//oo are determined using the RR. The division parameter value s = s* = 3.3783 
is defined by condition 

/i; = V^(r* + g)(M*) -1/2 = 0. 

The fixed point corresponds to the presence of the long-range interactions with the 
infinite correlation length. For describing the system behaviour, we should perform 
the infinite number of iterations (n 00). Nevertheless, the presence of the field 
or the temperature value distinguishing from the critical one, leads to the deviation 
from the fixed point {w*,r*,u*). In this case, the correlations become finite. At 
small values of the field h {h = ^ 0.01) and reduced temperature r = (T — Tc)/Tc 
(r = 0-i-O.Ol), i.e., at small deviations from the fixed point, we can use the linearized 
RR whose solutions are given by 

d/2. 



UJn = UJ* - So'^h'E^, 



4 ) 



Un = u* + ctl{f5^{Q,)fTi%l'\El^ + Ck2E2. (14) 

Here Ei = 20.977, E2 = 7.374, £'4 = 0.397 are the eigenvalues of the RG linear 
transformation matrix, T^f = u*^''^R2i{Ei - i?22)"\ T§ = M*"^/^i?42(^2 - Ru)~^- 
The coefficients c^^^ and 0^2 are expressed by relations 



$ - g - tS^uo^o '/'(/3$(0)/3e$(0))- 



-'24 



X 



7^(0)^(0) 
-^24 -^42 



Ck2 



no-n*-Ti°Vo^'/5<^>(0)(ro-r*: 



1 7^(0)^(0) 
i -^24 -^42 



where y^o = (g(l — s*~^)//oo)^. These quantities do not depend on the field. Using 
the condition of the hnearity of deviations from the fixed point, we can define the 
number of iterations n = Hp. This quantity determines the size of the scahng region. 
Taking into account the values of the eigenvalues Ei and coefficients in terms with 
the variables h and r in Eqs. fll4|l . one can see that the deviations from the fixed 
point are formed mainly in the first two expressions for a;„ and r„. Thus, the quantity 



Up is found using the equation 



-,np+lN 2 



(0). 



.*2 



(15) 



The right-hand side of Eq. ()15|) determines the maximal deviation, which satisfies 
the linearity condition. The choice of this quantity is also based on the condition 
related to the coefficient dn+i{k). When the block lattice size riches the correlation 
length, the coefficient dnj,+i{k) in expression (fTTj) in the case oiT > T^. changes sign 
and takes on the positive value. Therefore, it is possible to perform the next step of 
integration of the partition function in the Gaussian approximation. 

The quantity Up = Upij, h) (or the point of exit of the system from the scaling 
region) is a function of field and temperature. This characteristic provides the 
crossover between critical behaviour controlled by temperature and one defined by 
the field variable. The equality of the left-hand side terms in Eq. (fT3|) corresponds 
to the case described in the field-temperature plane by relation 



r oc 



(16) 



In the vicinity of values of the field and temperature satisfying relation ()16|). the 
scaling variable is of the order of unity and the series expansions in the scaling 
variable are not effective. It should be noted, that the RR (fT^ and their solutions 
(jl4j) are valid in the high-temperature region as well as in the case of T < at 



n < Hp. 



3 Critical exponents and thermodynamic charac- 
teristics of the system 

3.1 The case oi T > 

As was mentioned above, the block lattice size after the Hp iterations is commensu- 
rate with the correlation length. Increasing the number of iterations n and excluding 
"non-essential" variables pk, we renormalize the coefficients a\'^\ At n > Hp, due to 
changing the sign of the coefficient in the quadratic term in the exponent of relation 
(fTTj). we can calculate the quantity In+i appearing in ® on the basis of the Gaussian 



approximation. In the case of T > T^, the coefficient dn +i increases gradually and, 
for the quantitative results, one should perform the additional step of integration 
)12j . Then the partition function (jH)) takes on the form 

Z OC QoQl--Qnp+lInp+2, (17) 

where the quantity Inp+2 has the form similar to expression After perfor- 

ming this iteration, the quadratic term in Inp+2 dominates over other terms for 
all 7^ 0. Calculating the integral J„ +2 with respect to variables pk with indices 
< k < Bnp+2, we obtain 



n 



AO) 

'np+2- 



The quantity 



AO) 

'np+2 



(dpo) exp 



N'/'h'po - -dnp+2i0)pl 



(19) 



is crucial for system description by two reasons. Firstly, as a result of the RG 
transformations, every next block lattice includes the previous one and exhibits its 
physical properties. Since, the integration with respect to the variable po is the last 
step of calculation, this contribution is the most essential for the system description. 
The case of = corresponds to the total renormalization of the coefficients of the 
distribution of the order parameter fluctuations. Secondly, the average value of 
the variable po is the order parameter of the system [3 CHI- Thus, performing the 
substitution of the variable po = \/N po and evaluating the integral (fT^ by the 
steepest-descent method, we obtain the expression for the system free energy 

1 



k=0 



1 



h'a^s 



ip+2S 



-2^2 
(7n 



4! 



The quantity cxq is the root of the cubic equation 

_o 1 ^ . 



'r2S (To 



6 



^Mnp+2(To = 0, 



(20) 



(21) 



which is obtained using the extremum condition for the expression in the exponent 
of the relation (fT^ with respect to the variable po- Here, we use the substitution of 
variable 

(22) 



Po 



The choice of the root of Eq. (PT|) is based on the condition of the minimum for 
the free energy. The contributions related to variables pk with nonzero values of the 
wave vector are small in comparison with the contribution (|2(jp. The integration with 
respect to these variables is necessary for obtaining the renormalization coefficients 
appearing in expression ()2U|) . which is the system free energy at T > Tc. The relation 
corresponds to the equation of state. 



3.2 The case of T < 



Performing the integration of expression (fTT|) in the case of the low-temperature 
region, the displacement procedure proposed by Yuhnovskii (TU] is used. It is pro- 
vided by the substitution of the variables 



Pk = Vk + (ThVN6k. 



(23) 



After performing such a substitution, the coefficient in the quadratic term in Eq. 
(fTTjl becomes positive and dominating, as in the case of T > Tc. It takes on the 
form 



dh{k) = dnp+i{k) + 



(np+l) 2 



N 



(24) 



Up + l 



The detail procedure for obtaining this expression is given in [T^. Calculating the 
integral lup+i with respect to the variables ?7k at /c 7^ in the Gaussian approximation 
and returning to the variable Pq, we obtain 



n 



TT 



AO) 



(25) 



where 



AO) 

'np+l 



(dpo) exp 



(26) 



Performing the integration in expression for Inp+i by the same method as in the case 
of T > Tc, we find the system free energy 



k=0 



1 



and equation of state 

in the low-temperature region. 



CThS ~ 



(27) 
(28) 



3.3 Critical exponents and scaling functions 

The free energies (1201), (E3) equations of state (j^ . in the cases of T > T^. 
and T < Tc are explicit analytic expressions with coefficients, which depend on the 
microscopic parameters, in particular, on the lattice constant and parameters of the 
interaction potential. At h = 0, taking into account Eqs. ()15|) and (j^Hj) . we get the 
scaling representation 

pf,-) = B't^, T < T„ (29) 



where (3 = i^/2, B' = ah{c^ki /q)^ , and the critical exponent of the correlation length 
u is calculated by the formula 

In s , , 

0.609. (30) 



\nE< 



The value of the critical exponent u = 0.609 corresponds to the model app- 
roximation with 7] = 0. In pni; the quantity u = 0.643 is found for the similar 
approximation (the quartic potential but t] 0). As was shown in ^3^1^], the 
better convergent results within the approximation rj = can be obtained by using 
more complicated models. In particular, u = 0.637, u = 0.646 and u = 0.649 for 
models p^, and p^°, respectively The critical exponent of the correlation 

length for 0(1) model calculated using Wegner- Houghton RG equation takes on the 
values 0.527, 0.586, 0.673 and 0.759 for the second, third, fourth and fifth order field 
truncations (corresponding to 0^, 0^, 0^ and 0^°), respectively 
On the critical line r = at /i > 0, we have 

-pi^^ = D;y'h'K (31) 

where Dc = r*/ {cqSq^"^), and the critical exponent of the order parameter 5 is defined 
as 

*=iMl = 5. (32) 
In s 

Note that quantities ctq in Eq. ()22|) and ah in Eq. (PHj) coincide on the critical line. 
In order to compare the results with others, we use the explicit form for the scaling 
functions. In Fig. [H the equations of state and (PHj) are presented in the scaling 
form 

fciz) = Mh-^" (33) 

(the solid line), where z = i/h^^^^ is the scaling variable. Here h = h' /Dc, M = po, 
and t = tB'^/^ . Our results accord qualitatively with the well-known parametric 
representation of the equation of state 21j and Monte-Carlo (MC) results [22] (the 
dashed line). It should be noted that the results obtained in this work do not 
require any parametrization and allow one to take into account the influence of the 
microscopic parameters of the Hamiltonian (^on the critical amplitudes. 
In Fig. 121 the system susceptibility is given in the scaling form 

f^{z) = xDch}-'", (34) 

where x = dM/dh' . 

Some disagreement of our data with ones obtained in [21] can be connected 
with the model approximation, since present calculations are carried out using the 
simplest non-Gaussian model p'^. 



4 Conclusions 



The description of the ?>D Ising-like magnet behaviour in the vicinity of the critical 
point is proposed. The calculations are carried out by the CV method with using the 



RG transformations of the Wilson type. It is established that contributions related 
to nonzero values of the wave vector are small in comparison with one from the vari- 
able po; the average value of which is the system order parameter. The dependence 
of the limit magnitude of the long-range interactions on the temperature and field 
variables is presented. The scaling region size is determined as function of temper- 
ature and field by using Eq. (jl5|) . The universal and nonuniversal characteristics 
of the system are obtained simultaneously. Particularly, the analytic expressions 
for the free energy and equation of state, which are valid in the critical region in 
the whole field-temperature plane, are presented. The system susceptibility is con- 
structed in the scaling form. As is seen from Fig. |2l the susceptibility maximum 
location on the temperature scale is in good qualitative agreement with one obtained 
using the parametric representation of the equation of state and MC simulations. 

The calculations of system nonuniversal characteristics are based on the physi- 
cally well-grounded and mathematically rigorous use of the phase space of CV. They 
are performed starting from the Hamiltonian, which contains the microscopic para- 
meters. The results are obtained using the quartic approximation for the effective 
Hamiltonian [p^ model) without taking into account the correction for the potential 
averaging (analog of the LPA). Due to removing the cubic term by making the shift 
of CV, the presence of the field is represented only by the linear term. This result 
is valid also for models p^, p^, etc. In other words, the presence of the field caused 
to appearance of the additional relation for Un in Eq. (fT^ . 

The similar estimations can be generalized to the models with an n-component 
order parameter 
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Figure 1: The scaling function of the order parameter. The solid line corresponds to 
Eqs. (1221) at T > Tc and T < Tc, respectively. The dashed line corresponds 

to the parametric representation of the equation of state |2I] and MC simulations 




Figure 2: The scaling function of the susceptibihty (the sohd hne is the CV method 
result, the dashed line is the result obtained in [211 1221 )■ 



